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The correlations of the electrons field-emitted from a superconductor are fully analyzed, both in 
space and time. It is proposed that a coincidence experiment would reveal a positive correlation 
between the electrons emitted in opposite directions. The electrons can be entangled and can even 
violate Bell's inequality. The crucial role played by Andreev's process is scrutinized, analytical 
formulas are derived for the correlations, and the physics behind the phenomenon is clarified. 
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I. INTRODUCTION 

Coherence is a fundamental notion in quantum me- 
chanics. It is related to the superposition principle and 
plays a fundamental role whenever intrinsically quantum 
phenomena take place. When many identical particles 
are involved, the global picture can become very rich: 
quantum coherence, together with quantum statistics, 
may drive the system into highly nontrivial states, en- 
dowed with a variety of interesting features. The su- 
perconducting state of electrons in solids is one of such 
examples. It exhibits very interesting physics and has at- 
tracted the attention of many researchers from different 
perspectivesji 

The superconducting state is fully characterized by a 
set of correlation functions. If some particles are emitted 
from the superconductor, they retain some features of 
the correlations in the source. A superconductor is char- 
acterized by its long-range coherence, both in space and 
time, and by Cooper-pair correlations. The long coher- 
ence can yield a well-monochromatizcd beam of electrons, 
which can greatly improve the quality of a microscopci^ 
When a Cooper pair is emitted, its singlet spin state 
can be a useful resource of entanglement, with applica- 
tions in quantum informatioui^i^i^i^iiiSiiSiiiS It is also in- 
teresting to regard the emitted particles as probes of the 
sourceJ^ Their correlation functions reflect the features 
of the source and can behave in remarkable ways. Posi- 
tive correlations among electrons have been discussed in 
the context of superconductors 

A coincidence experiment is a very direct way to detect 
two-particle correlations. Recent technological progress 
has made possible coincidence experiments in a variety of 
systems li^ii^ii^ii^iii Among these experiments, we shall 
focus on i) a coincidence experiment in field emissionj^^ 
ii) the observation^ of the field emission spectrum from 
a superconductor.— A combination of these two experi- 
ments can lead to a challenging possibility: a coincidence 
experiment in the field emission from a superconductor. 

The nonlocality of the electrons field-emitted from a 
superconductor has been recently analyzed.— This arti- 
cle provides a thorough explanation of this phenomenon, 
scrutinizing several additional effects and the role played 



by some genuine superconductivity phenomena, such as 
the Andreev's process. The emission process will be dy- 
namically described in the framework of quantum field 
theory and the beam profile will be naturally prepared 
by the dynamics itselfj^ It is crucial to work in 3D space 
to capture the Cooper-pair correlations. This also en- 
ables us to discuss the lateral coherence length<^ Con- 
cise and useful analytical formulas will be derived, which 
will clarify several facets of the physics behind the phe- 
nomenon. The present analysis also includes correla- 
tions in time. The Andreev process will be shown to 
play a crucial role for these peculiar correlations and 
nonlocality^'^'^'i'^^^ 

The article is organized as follows. The problem is set 
up and the Hamiltonian of the system is given in Sec. [TTl 
The dynamics of the emission is solved and the stationary 
beam of electrons is obtained in Sec. IIIII The correlation 
functions and their spectra are computed in the station- 
ary beam and some analytical formulas, whose structures 
clarify the underlying physical processes, are presented in 
Sec. IIVI The coincidence spectrum is analyzed in detail 
in Sec. |Vl where the effects of the superconductivity of 
the emitter arc disclosed. Section I VII is devoted to the 
entanglement and nonlocality of the emitted electrons. 
The robustness of the correlations due to the supercon- 
ductivity of the emitter is discussed in Sec. IVIII Finally, 
the whole analysis is summarized in Sec. I Villi and some 
details are presented in Appendices EHDl 



II. SETUP 

A. Hamiltonian 

Let us start off by discussing electron emission from a 
superconductor into vacuum. We focus on the emission 
from a nanotip, see Fig. [TJ To this end, we set up the 
following Hamiltoniani^ii2ii^i2i in 3D space: 

H = Ho + XHt, Ho=Hs + Hv, (2.1) 
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FIG. 1: Emission of electrons from a superconductor in 3D 
space and detection of two electrons emitted in different di- 
rections. 



where 



M^0j(r)^|(r)0t(r-)0i(r) 



(2.2) 



is the Hamiltonian of the superconducting emitter™ and 
Hv^Y. I d'ri^lir) (-^V^) ^,(r) (2.3) 

is that of the electrons propagating in vacuum, s denotes 
spin and we set h= 1. (j)s{r) and V's(^) a-i'c the fermionic 
field operators satisfying the canonical anticommutation 
relations 



{Ur),4>l,{r')}^{^s{r).i^l{r')}^5,s'5^{r-r') 

(2.4) 

with other anticommutators vanishing. The Coulomb re- 
pulsion is neglected. 

Ht describes electron tunneling through a potential 
barrier surrounding the emitting region, and A charac- 
terizes the strength of the tunneling transmission. Ht is 
given byi^^>20^ 

Ht=Y1 [ [ (Tpkcl.aus + T;kal,c 

where aks and Cp^ are the annihilation operators in mo- 
mentum space of the electrons inside and outside the 
emitter, respectively, and are related to the fields in con- 
figuration space by 



Mr) 



CpgC 



ipr 



(2.6a) 
(2.6b) 



Tpfc are the tunneling matrix elements, for which we 
takoiS^ 

Tpk = {p\h{p)g{r)\k) 



Hp) 



(fir 
(27r)3 



g{r)e 



^h[p)~g{p-k). 

(2.7) 



That is, an electron with momentum k in the emitter is 
annihilated by aks^ filtered by g{r) and h{p), and emit- 
ted outside with momentum p by c)^^. The function g[r) 
specifies the emitting region and |/i(p)P represents the 
momentum (energy) dependence of the tunneling proba- 
bility through the potential barrier surrounding the emit- 
ting region. In this article, we consider a simple spheri- 
cally symmetric setup with 



^(27ru;2)3 



m 



1 



(27r)3 



h{p) = 



(2.8a) 
(2.8b) 



where w characterizes the size of the emitting region, Sp is 
the energy of an electron, to be given below in (|2.18p . and 
Ec controls the low-energy cutoff of the tunneling spec- 
trum. The function h{p) chosen here might not always 
accurately describe the tunneling probability through the 
potential barrier (for which the relevant quantity might 
be the energy related to the motion normal to the poten- 
tial surface, instead of the total energy Sp). In the follow- 
ing discussion, however, we are interested in the far field, 
for which the above choice will turn out to be appropri- 
ate, since only the momentum normal to the surface is 
relevant to the far field. 



B. Mean-Field Approximation 

Let us introduce 



where 



no = Ho- ijN, 



N^Ns + Ni, 



(2.9a) 
(2.9b) 

(2.10) 



is the sum of the numbers of electrons inside (Ns) and 
outside (Ny) the emitter, defined respectively by 



Ns - 



^ / d'r4{r)Mr), (2.11a) 
J2 J d'r^pl{r)Mr)- (2.11b) 



Notice that N is a constant of motion, i.e. [H,N] = 0. 
The time-evolution operator can therefore be split as 



iHt 



i/iAftg int^ g^j-^j exponential factors like e" 



ifj,t 



factorize away from the quantities of interest. Indeed, 
the Heisenberg operators are factorized as 

Vs(r,t) = e'^Vs(r)e^'^* = ^,(r, ^e"'''*, (2.12) 
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and so are the correlation functions, 



('0li(n,tl)V's2(''2,t2)) 



(^I^(ri,ii)^,,(r2,i2)>e'''(*i-*^), etc., (2.13) 



where 



iHt 



-mt 



(2.14) 



describes the dynamics of the field in the picture intro- 
duced by the unitary transformation e*'^^*. 

We shall work in such a picture, with ^ the Fermi level 
of the superconducting emitter. This choice is conve- 
nient for the mean-field approximation, which enables 
one to diagonalize the Hamiltonian Tig via the Bogoli- 
ubov transformation^ 



(2.15) 




-ki/ 



with 



Uk 



V2 



„iS 



£fc_ 



V2V ^fe' 



(2.16) 



to get 

Hs = Hs- ^lNs = 

Hv =Hv - uNv 

where 

fi, 



J2 / d^kiJkal^aks 

s=hl 

J d^pEpcl.Cps, 



(2.17a) 
(2.17b) 



" 2m 



^k 



are the energies of an emitted electron in vacuum and of a 
quasiparticle excitation in the superconducting emitter, 
respectively, measured relative to the Fermi level of the 
emitter, and 



(2.19) 



is the gap parameter of the superconductor. Throughout 
this article, kp = y/2m^ and Xp = 2-K/kF are the Fermi 
momentum and the Fermi wavelength, respectively. 

The state of the superconductor is characterized by the 
Fermi distribution of the quasiparticle excitations. 



f{ujk)5ss'5''{k~k') 



with 



1 



(2.20) 



(2.21) 



where T is the temperature of the emitter and fc^ the 
Boltzmann constant. The spatial extension of a Cooper 
pair is characterized by the correlation length of the two- 
point correlation function ((/)|(ri)(/)j^(r2)) in the super- 
conductor, which is given (at zero temperature) by Pip- 
pard's length^ 



kp 



2fi 



TrmlAI TrfcFlAI 



C. Coincidence 



(2.22) 



At the initial time t ~ 0, the emitter is in the supercon- 
ducting state at a given temperature T and the outside 
of the emitter is vacuum. Such an initial state is a prod- 
uct state of the superconducting state characterized by 
(|2.20p and the vacuum state for Cps- Starting from this 
initial condition, the electrons start to tunnel out of the 
emitter into vacuum. The emission process is dynami- 
cally described according to the Hamiltonian (|2.9p with 
p.5|) and ()2.17|) in the mean-field approximation. Af- 
ter a certain transient period, the emission approaches a 
nonequilibrium steady state (NESS)i^ 

We shall count coincidences, in the NESS prepared in 
this way, between two detectors located at 



ri,2 = (±rsin(0/2),O,rcos(0/2)), 



(2.23) 



at the same distance from the emitter but in different 
directions, as depicted in Fig. [TJ The numbers of counts 
of one- and two-particle detections, irrespectively of the 
spin state, are proportional to 



p(i)(r,t) 



J2 (V'l(r,t)V.(r,i)), 



(2.24) 



(2.18) p'^^\r2,t2;ri,ti) 



(i2 >ii), 
(2.25) 

respectively, where the average is taken over the initial 
state. The deviation of the normalized two-particle dis- 
tribution function 

^/ , , N p''^Hr2,t2;ri,ti) /, ^ , n 

Q{r2, t2] ri, ti = -J-- — 7YT7 ^2 > h) 

pW{r2,t2)p^^>{ri,ti) 

(2.26) 

from unity reveals correlations between the two electrons. 



III. DYNAMICS OF EMISSION 



Let us solve the dynamics of the emission. The Heisen- 



berg equations of motion for Cps (t) 



Cr\aC 



-mt 



etc., 
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read 



where Qppi (t) is given by the inverse Laplace transform 



dt i?t 



'-Pi/ 



"pi. 



iX / d-'kTr 



a' 



(3.1a) 



where 



\ 




° 















TpkUk 
"-^p(-fc)^'fc 



-T, 



p(-k)Vk 

'^pk'^k / 



This set of equations is solved by means of the Laplace 
transform to yield 



Spi< 



Cp'T 



dV^pp'Wi ; 

di! ^ d^k[g{t~t')T] 



, I oik] 



(3.4) 



Qpp'it)^ ^^Gpp-{s)e^' 



(3.5) 



(3.1b) of 

(3.2) 
(3.3) 



Qpp'i-'^) = (-^ + ^£v)^\p - P') + X^ICpp'is), 



ICpp'{s) ^ / d^kTpk 



(3.6) 
(3.7) 



with Cb the Bromwich path running parallel at the 
right of the imaginary axis of s, and [Q{t — t')T]pk = 
J d^p' Gpp'{t — t')Tpik like a matrix product. 

In order to obtain the NESS, it is convenient to 
move to the interaction picture defined by Cps(t) = 
e^'^°*°Cps(t)e'^''*° , since the initial state is invariant 
under this transformation. In this picture, Qppi (t) = 
Qpp'{t)e^^p'^° asymptotically behaves, for i,to oo, 
keeping t — to finite, as (Appendix IX]) 



Gpp'it) ^ e-'^'-^^-^'^S^p ^ p') - X 



' A'",(— te„,+0+) 
2 I i{ep-V)+0+ 



where 



^pp' (^) 



pp' ^ P 



-j(ep+£p,)+0+ -*(<rp-S')+0+'' 



+ 0(A4), (3., 



^Pk^p'k [s+iuk ^ s-it^fc 



^^pfc^p'(-fc)"fe"fe (s-Lfc s+Lfc) 



TpkTpH_k)UkVk 



\ s—iu. 



(3.9) 

Note that we need to retain up to O(A^) to collect the lowest-order contributions to the correlation fimction (|2.25p . 
which is 0{X'^). Up to this order, the Heiscnberg operators p.4p in the NESS in the interaction picture read 



d'p'gpp>it)[, )-x d^ki 



/'^fcT\ 

e-'^^-i^-^o) +0{X^). (3.10) 

^_ , ^ , \"L/ 



'p'i ' ' 



IV. CORRELATION FUNCTIONS 
A. Structures of the Correlation Functions 



The correlation functions of the emitted electrons in the NESS are now readily computed. Let us start with 
two-point correlation functions. In the limit ti, t2 oo, keeping ti — t2 finite, one gets, up to O(A^), 

l{ri,ti;r2,t2) = (■(/i|(ri, ii)-0|(r2, i2)) = (■^{(ri, ii)'0x(r2, i2)) 
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(4.1) 



where 



Xo{ri,ti;r2,t2) + Xthiri,ti;r2,t2), 



Xth{ri,ti;r2,t2) = -A^ ^ / d^fe (ri, ti).^fc^(r2, t2), 



(4.2) 



(4.3) 



Xo{ri,ti;r2,t2) 
= AM d^kukVk 



y / d^kukVk 



d^pi d^p2 Tp,kTp,^^k)e'P''-'e'P'---^ 
^/i2nf ^/i2nf (cpi - i^k - «0+)(ep, + Wfe - iO+) 

I 



JPl -ri jp2 •I-2 p-lEpi (ti -(2) 



(4.4) 



while all other two-point correlation functions vanish. In 
these expressions, 



Vk-{r,t) 



d^p 



TpkUk 



^(27r)3 i Sp - ujk 



d^p 



TpkVk 



^i(p-r-ujkt) 

(4.5a) 



y^(2nf i Sp + LUk - iO+ 



(4.5b) 



are the wave functions of the emitted electrons in the vac- 
uum originating from the quasiparticle excitations above 
and below the Fermi level, respectively, and they propa- 
gate asymptotically as spherical waves (Appendix [B| 



ipk+ir,t) ~ mV^jrT, 



1 



ipk-ir,t) ~ mv2'K9{fj. - uJk)T(^p^r)kVk—e' 



for kpr ^ 1, where 

p± = p{±iUk), p{E) = V2m(M + E) 



(4.6a) 
(4.6b) 
(4.7) 



are the momenta of the emitted electrons correspond- 
ing to the energies ±0;^, and 6{x) is the Heaviside theta 
function. 

The correlation function 7 in (j4.ip is a mixture of 
the wave functions </3fc±(T,t) over the Fermi distribu- 
tion /(zbtjfc) in the emitter, with the relevant densities 
of states of the quasiparticle excitations, and 
The quasiparticles are emitted outside and the electrons 
propagate with the wave functions 'fk±{f,t)- The cor- 
relation function 7 is essentially a one-particle density 
matrix and describes the one-particle state of the emit- 
ted electrons. 



The other correlation function x in (|4.2p - (|4.4p . on the 
other hand, describes pair emission. It contains the prod- 
uct UkVk (= A/2u;fe) and is responsible for the effects of 
the Cooper-pair correlation in the superconducting emit- 
ter. Indeed, xo describes the emission of a couple at zero 
temperature, from the Cooper pairs, and takes into ac- 
count the Andreev process . '^'^'^'^ Notice that the last con- 
tribution in (14. 4p vanishes on the energy shell Sp^ = uJk- 
The main contribution is the first addendum and it is 
exactly the amplitude for the emission of a couple [see 
(|C7P in Appendix [C] . The contribution of the thermal 
excitations at finite temperature is taken into account by 
Xth, which vanishes at zero temperature. 

Since the solution of the Heisenberg equations of mo- 
tion, (|3.4p . is linear in the fermionic operators, the one- 
and two-particle distribution functions, (|2.24p and (|2.25p 
are both given in terms of the two-point correlation func- 
tions: the former is given by 



pW(^,i) = 27(r, i;r,t). 



(4.8) 



while the latter is cast through Wick's theorem into 

P^'^(r2,i2;n,ti) 

= 47(r2,t2;r2,t2)7(''i7ii;''i7*i) 

- 2|7(r2, t2; riM)? + 2|x(r2, i2; ViM)? . (4.9) 



The normalized two-particle distribution ()2.26p is there- 
fore given by 



g(2;l) = l 



l7(2;l)P 



lx(2;i)|^ 



27(2; 2)7(1; 1) 27(2; 2)7(1; 1) 



, (4.10) 



where the arguments (r2, t2; Ti, ti) are abbreviated to 

(2;i). 

The second term of the normalized coincidence (j4.10p 
is responsible for the antibunching of electrons (reduction 
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of the coincidence rate when the two detectors are close 
together and the detection delay time is small), while the 
third gives a positive contribution to the coincidences. It 
is instructive to rewrite the two-particle distribution (|4.9p 



J 



in terms of the two-particle wave functions: by substitut- 
ing the two-point correlation function 7 given in (|4.ip . 
i.e. the contributions of the quasiparticle excitations, the 
first two terms in (|4.9p read^S 



(72 .CTi— ± 



<^..fc,.,(2;l)'+ <^,,fc,.,(2;l) V2|x(2;l)|^ (4.11) 



r 



where 



(± 



fclCl.fc20-2(-'-' 

= {vk^a^ {^)'^k^<y^ (2) ± {^)'fk^a^ (2)) 

(4.12) 



are the symmetrized/antisymmetrized two-particle wave 
functions of the emitted electrons. This expression clari- 
fies that 3/4 of the contributions of the quasiparticle exci- 
tations are provided by the antisymmetric wave function 
while the remaining 1/4 is given by the symmet- 
ric one Recall that the state of the fermions as 
a whole should be antisymmetric. The symmetric wave 
function in space therefore corresponds to the antisym- 
metric state in spin, and vice versa. As already men- 
tioned, 7 describes the emission of the single particles. 
The electrons in the emitter are not spin polarized when 
one looks at the single electrons and they appear to be 
equally emitted from the triplet (antisymmetric in space) 
and singlet (symmetric in space) spin states. This just 
gives the background: the interplay between these contri- 
butions yields antibunchingi^'22 but does not reveal the 
effect of the Cooper-pair correlation in the emitter. 

The Cooper pairs do not significantly contribute to an- 
tibimching, but yield the positive correlation, i.e. the last 
contribution in (|4.10p . Let us look at the Andreev ampli- 
tude, the first contribution to xo in (|4.4p . This expression 
entails that one electron in a pair originates with momen- 
tum k in the emitter, while the other with momentum 
—k. Notice here that the integrand is symmetric under 
exchange between k and — fc and the two alternatives are 
symmetrically superposed. This symmetric superposi- 
tion reflects the fact that the Cooper pairs in the emitter 
form singlet couples, antisymmetric in spin. This symme- 
try in the amplitude with respect to the spatial degrees 
of freedom gives rise to bunching: two electrons bunch 
when they originate from momenta with the same mag- 
nitudes but oriented in opposite directions. This is the 
origin of the positive contribution in (j4.10p . 

Positive correlation of electrons due to superconduc- 
tivity has been discussed in the literature, in the cur- 
rent fluctuations in transport setupsJ^ For instance, it 
is shown that the cross-correlation of the currents at two 
different normal leads connected to a superconductor can 



be both negative and positive, while it is quite generally 
negative in normal devices. It is pointed out that the pos- 
itive correlation is due to the processes involving "crossed 
Andreev reflections," which provide transport of particles 
with opposite charges from one normal lead to the other 
via the superconductor (an electron in one normal lead is 
converted into a hole in the other normal lead, and vice 
versa) ^1516,7,8 ^10,12^24,25 ^^y^ gjgj^ q£ ^^^^ overall cross- 
correlation is determined by the competition between the 
negative contributions due to the normal scattering and 
the positive contributions due to the Andreev reflections. 
In the present field-emission setup, there is no "hole" in 
vacuum, but the crossed Andreev reflection is equivalent 
to the emission of electrons from a common Cooper pair 
toward different detectors. In the present analysis, it is 
clear from the formula, as just explained above, that the 
positive correlation is due to bunching as a result of the 
symmetric wave function of a couple of electrons emitted 
from a Cooper pair, reflecting its singlet spin state. In 
addition, it is shown below that the negative and posi- 
tive contributions appear separately in the present field- 
emission setup, and the pure positive correlation is ob- 
served: one need not argue the competition between the 
negative and positive contributions. 

The amplitude for the Andreev emission (|4.4p is also 
symmetric under the exchange between ri and r2. This 
symmetry reveals that the emitted electron pairs are also 
in the singlet state. This point is explicit in the expres- 
sion (|C7p . The entanglement of the emitted pair is one 
of the main subjects of this article and is explored below. 

It is interesting to note that such emission of the singlet 
pair is not simply the direct emission of the Cooper pair 
since virtual processes are involved, as is clear from the 
expression in (jC7[) . 



B. Spectral Representations 

Let us look at the energy spectra of the cor- 
relation functions, r(ri,r2;-E), A'th (t^i , T2 ; ^') , and 
'^o{i'i,f2]E), which are related to the correla- 
tion functions 7(''i: *i; ''2, ^2), Xth(T'i, ii; r2, ^2), and 
Xo{i'i-,ti\r2,t2) in (|4.ip - (|4.4p . respectively, through the 
Fourier transformations, 
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r dE , 
7(ri,ti;r2,t2) = / —T{ri,r2]E)e 



iE{ti-t2) 



(4.13) 



etc. By plugging the asymptotic forms of the wave functions in ()4.6p together with the tunnehng matrix ()2.7p - (j2.8p . 
the energy spectra in the far field region {kpr > 1) in the NESS {t oo) read (Appendix [B| 



T{r,,r2;E)^27rAe{\E\-\A\) 



Xa-,{ri,r2;E)^27rA0{\E\-\A\) 



\E\ 



A 



Z{E)f{E)e^^^<^, 



Xoiri,r2;E) ~ ttA 
for \E\ < n, where A = X"^ m"^ /2{2t:)^ w'^r'^ , and 
ZiE)^ 



A 



Zth{E)f{\E\)e'^P^'^'>+P^-^'>^'' 



fZoiE)e 



i\p(E)+p{-E)]r 



(4.14a) 
(4.14b) 
(4.14c) 



2cos(0/2) V ^ / ^ 



- 2g(i?) ' 

Zo{E)=Ztb{E) + 5Zo{E), 



(E)] 



(T = ± 



-«.'[fc<,(i5)-g(i5)]' _ -w^lk„{E)+q(E)]'' 



SZoiE) = I 



2g(i?) 



5] a (siwiKiE) ~ q{E)]) - ^(wiKiE) + qiE)]) 



(4.15a) 
(4.15b) 

(4.15c) 



r 



with p{E) defined in (|4.7p and 



of the spectrum is controlled by the factor f{E)e^^^'^ 
and the spectrum spans the range 



k±{E) = J2m ^ ± v/-E2 - |A 



q{E) = ^[kl -p{E)p{-E) cos0]/2. 
Note that 

S(z) = e"""^ erfi(z) 



(4.16) 
(4.17) 

(4.18) 



- Ec <E <il/kBT-l/Ec)-\ 



(4.20) 



A remarkable feature of the spectrum P{E) shown in 
Fig. [2] is the presence of the gap for the superconducting 
emitter: no emission from the quasiparticles in the range 
— |A| < E < |A|. In addition, the spectrum diverges at 



is a slowly varying function of z and t hat we h ave 
chosen the branch of the square root ^/E"^ — | Ap = 
V|AP - -£2 for \E\ < |A|. 

The energy spectrum of the emitted electrons, 



P{E) (X ATrr'^T{r,r;E), 



(4.19) " 



ST 



|A|= 0.003 M 
^ 0.001 
^ 0.000 

S-H 



kBT= 0.2 [10- V' 
Ec = 4.5 [10- V' 
w = 5[Ai?] 




is shown in Fig. [21 where on the basis of the formulas 
(|4.14ap and (|4.15a|) . The low-energy cutoff reflects the 
energy dependence of the tunneling probability: electrons 
with a lower energy feel a higher and thicker potential 
barrier and tunneling is suppressed. In the present model 
with (|2.8p . the low-energy slope is characterized by the 
factor e^Z-^c, The high -energy cutoff, on the other hand, 
is due to the Fermi distribution f{E) and is therefore FIG. 2: (Color onl ine) E nergy spectr um of the emitted elec- 
mainly controlled by the temperature. The band width ^=^0^^' M ^i^h gUil) and glSSJ. 
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FIG. 3: (Color online) Spectrum of the Andreev emission for 
e = TV with (|4.14b|l - (|4.14c|l and (|4.15b|) - (|4.15cp . The param- 
eters are the same as those in Fig. [2] 



the edges of the gap, reflecting the density of states of 
the quasiparticles in the superconductor. 

It should be noted, however, that we are looking at 
the two-point correlation functions at the lowest order 
in the tunneling Hamiltonian, i.e. up to O(A^) [see (j4.ip - 
(|4.4p ]. The correlation function 7 in (|4.1|) up to this order 
counts only the direct emission from the quasiparticlc ex- 
citations, but there also exist emissions from the Cooper 
pairs. Let us look at the spectrum of the other two- 
point correlation function, which describes the An- 
dreev emission from the Cooper pairs. Figure[3]shows the 



spectrum for ri = — r2, i.e., when the pair electrons are 
detected in opposite directions. It exhibits nonzero spec- 
trum even inside the gap of the quasiparticlc spectrum. 
This is the distinctive feature of the Andreev process. 
If one collects higher-order processes, the single-particle 
spectrum P{E) also accounts for the contribution of the 
Andreev emission, and one sees a spectrum also in the 
gapi^ 

Suppose now that 



M > -Be (> fcisT", usually), 



(4.21) 



which is usually the case in actual experiments. Then, 
only energies close to the Fermi level \E\/ <C 1 coir- 
tribute to the correlation function 7. As for the other 
spectra and in (j4.14bp and (|4.14cp . the factor 



(4.22) 



rapidly oscillates for kpr ^ 1 and its stationary phase 
at the Fermi level E ~ Q provides the most significant 
contribution. The relevant energy range in the saddle- 
point approximation is 



\E\hL < V'^/kFT « 1. 



(4.23) 



Therefore, only energies close to the Fermi level \E\/fi <C 
1 are relevant to the correlation functions. In this regime, 
Z{E), ZthiE), and SZo{E) in (|4.15p are approximated by 



Z(E) ~ ^ ( ^~4:k:irw'^ s\n^ {e/4){l+E / fi~\A\^ /2fj.E) _ -ik^w'^ cos^ {6 /4){1+E/ fj,-\A\'^ /2fiE)\ 

^ ' cos(6'/2) V )'■ 

Zth(B^)e'[P(^'+P(~-^'l'' ~ /'e-4fcj,t«^ sin=[(u-e)/4]p(»V'r'?') sin(e/2)^-j£;V2K + 

sin(0/2) V 

<5Zo(-B)e'[P(^)+f(-^)l'' 



sin(6'/2) 



ikpwt^ y/E"^ - |A|2 e^^'^'^'' 
Jl sin(6'/2)' 

I 



(4.24a) 

(4.24b) 
(4.24c) 



where 



Kj^ = -j{r/2kF T ivw^), i/ = sin - - - cos - cot -, 

(4.25a) 



1 



Kf, ^ = —^(r/kp + iw^ cost 



(4.25b) 



Additional conditions have been assumed: 



l/2fci.rsin2(6i/2) < 1, ^kpw'^j^r^ 1 (4.26) 



for (|4.24bp and (|4.24cp . Note that Eq. (|4.24ap is still 
exact for a normal emitter A = 0. 



A 



^' Ikpw sin 



— ^' ( 2kpw COS' 



(4.25c) 



C. Analytical Formulas for the Correlation 
Functions 



Several concise and useful analytical formulas arc now 
available. In the following, let t = ti — t2- 
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1. Normal emitter A = at T = 2. Normal emitter A = at T > 

111 this case, x vanishes, while the inverse Fourier trans- 
form of (|4.14a|) with (|4.15ap [or equivalently, with (|4.24ap ] Even at finite temperature T > 0, if 
yields 



7(1; 2) 



cos{e/2) \Ec^ +iT~4{kj,wy^i)sm\d/4) l/keT > 1/ Ec - 4.kj,wy fi > and e"^/-^^ < 1 

g-4feJ.«,2cos^(e/4) _ g-p(£;-i+jr) \ (4.28) 

— ; — — — — ^ j . are satisfied, an analytical formula is available: the lower 

+iT - A{kpW I \i) cos (6'/4) / ^]^^ integration range in energy can be extended 

(^■2''') to —n — oo and one gets 



J 



cos(6i/2) \sm{'KkBT\Ec^ + ir - A{klw'^ / ii) sin2(6'/4)]} sm{-KkBT[Ec^ + ir - 4{kj,w^/ii) cos2(6'/4)]} J 

(4.29) 

I 

by making use of the formula (jDl|l . G{ri,r2; t) 

= y°° g Zth(i?)e'[P(^)+f(-^)l'-e'^^ 



g-4fc>=sin^[(77-e)/4]g(«,V'r « )sin(0/2) 



3. For superconducting emitter |A| > at T = 
The inverse Fourier transform of (|4.14ap with (|4.24ap sin(6'/2) 

2ni 

IAI?C, ClAirRTi 4-,V^^ 

cos(6i/2) 

^g-4fcp«)^sin^(0/4) _ g-4fc^u)^cos^(6l/4)^ 



yields / K+ 

7(1; 2) ^ --A_|A|Xi(|A|(ii;-i +*r)) (4.33) 

and 

X ( g-.-F«'^sin^(0/4) _ g-4fc^u)=cos2(e/-^^^ 

(4.30) 
<5xo(l;2) 

bv neglecting 0(En / fi) contributions, and that of ( |4.14cp . a 

with (I4.24bp and KMBi gives =-/ dE AZn( K)p^[p(^)+p(--^)1''p^-^^ 

V 27ri 

(4.34) 



dr' F{t - T')G(ri, rs; r') + 6xo{l; 2) ^ V^!^ Mi^g2.fcrrg-4«,^ co.^(e/2) /^g«orV2 

(4 3^^ sin(0/2) V2^z 

with 



-oo 27r -^£^2 _ I Ap 2i For the coincident detections with t = ti — t2 = 0, an- 

(4.32) other expression is convenient: 



J 



^0(1;2) ~ ^^^e"^"" L-44»^sin^[(.-e)/4]g-,r/2.^fc.?=g-(»V-^«^)[™(('/2)]^(2)((,^y,2 „ ^/^kF)/^^^) 



4i sin(6'/2) 



TT^y ir /kpw"^ 



(4.35) 



In (|4.33p and (|4.35p . the second contribution in (|4.24bp has been omitted by assuming e 2feJ-Mi^ ^ 1. As for the Bessel 
functions Ky{z) and H^\z), sec Appendix [Dl Note that A defined in (|4.25cp is almost constant A ~ 1 for w > Ai? 
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kBT= 

Ec = 4.5 [10- V] 
w = 5[Xp] 
r = 10* [Xp] 



|A| = 0.003 [H 
0.000 




r 



0' • • • ' • • • ' • • • ' • • • ' • • • ' 

0.0 0.2 0.4 0.6 0.8 1.0 

0/tt 

FIG. 4: (Color online) Normalized coincidence Q vs 6 for 
ti = t2, for normal and superconducting emitters. 



and 6 t:. 

An analytical formula is not available for the supercon- 
ducting emitter |A| > at a finite temperature T > 0. 
In the present analysis, the effects of temperature are 
taken into account only though the Fermi distribution 
function f{E), which controls the high-energy tails in 
the spectra. They are usually small, with Ec 3> ksT 
like in Fig. [2l even at room temperature. For instance, 
irkBT/ sm{TrkBT/Ec) ^ Ec in formula (|4?29l) (almost no 
temperature effect), and the other analytical formulas for 
T = agree well with the numerical estimations for the 
finite-temperature case with the value chosen in Fig. [21 
It therefore suffices in the following discussion, to use the 
formulas for T = 0. Physically, the gap parameter A is 
a function of temperature T, so that temperature effects 
appear through the gap parameter |A|. 



V. ANTIBUNCHING AND BUNCHING 

We are now ready to discuss the correlation between 
the electrons emitted from a superconductor. The nor- 
malized coincidence Q is shown in Fig.|3]as a function of 9 
at ti = t2 = t, for normal and superconducting emitters. 



A. Antibunching at 6 ~ 

The dip around ~ is due to the antibunching of the 
electrons: there are fewer chances to detect two electrons 
at the same time than expected from the counting rate 
of the independent detections of single electrons. The 
dip is Q ~ 0.5, irrespectively of the parameters. This 
is clear from formula (|4.10p for Q. Figure [3] shows that 
the superconductivity of the emitter does not affect the 
antibunching and x — for 6* ~ 0. Hence, formula (|4.10[) 
yields Q ~ 0.5 for 6* = and ti = t2, irrespectively 
of the details of the correlation function 7. Equation 
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y Jl 











FIG. 5: (Color online) Normalized coincidence Q vs 6 for 
ti = t2, for a superconducting emitter with different values of 
w. The other parameters are the same as those in Fig. |4]with 
lAI = 0.003 [ul. 



(|4.1ip shows that three antibunching contributions plus 
one bunching yields 1 - 3/4 + 1/4 = 0.5-ii^ 

The analytical formulas for 7 in Sec. IIVCI show 
that the width of the dip is governed by the factor 
g-8fcp u) sin (e/4) controlled by the size of the emit- 

ting region, w. The smaller the size of the emitting re- 
gion, the longer the lateral coherence length of the emit- 
ted electrons, and the wider the antibunching dipi^ See 
Fig. [SJ where the effect of w is scrutinized. 



B. Bunching at 9 ~ n 

The effect of the superconducting emitter manifests 
itself when the electrons are detected in opposite direc- 
tions, 9 ^ TT. Since a Cooper pair is formed by two elec- 
trons with opposite momenta k and — fc, they are emit- 
ted in opposite directions and exhibit strong correlation 
at 6 ^ TT. The spins of the Cooper pair are in the singlet 
state, and the corresponding symmetric wave function in 
momentum space gives rise to bunching, namely a posi- 
tive correlation, as already discussed in Sec. IIV Al 



The analytical formulas for x in Sec. IIV CI show that 
the width of the bunching peak is characterized by 
the factor e-sfeFiw-sm [(n-e)/4]^ ^^^^ ^^le lateral coherence 
length is the same as that for the antibunching. The 
Cooper pair does not affect the lateral coherence. 

The height of the bunching peak ^Qpcak = 
Q{r, t, — r, <) — 1 is, on the other hand, given by 



lx(r-,t;-r-,t)|- 
2|7(r,t;r,i)|2 
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?,2Kim/Ec) 



(5.1) 



the ratio of (|4.35p to (|4.30p . This is the ratio of the con- 
tribution of the pair emission to that of the single-particle 
emission, the background. The main contribution to the 
former is provided by the term represented by the Hankel 

(2) 

function when w ^ ^: it is greater than the other 
by a factor of order ^/w [see for instance the asymptotic 
form of the Hankel function in (jD9p ]. This pair emission 
is controlled by the parameters r/2kp^^ and w^/^^, both 
measured against Pippard's length ^, while the single- 
particle emission by Ec- In particular, the emission 
from a smaller emitting region results in a stronger cor- 
relation. See Fig. [3 which clarifies how the size of the 
emitting region, w, affects the bunching peak. 



C. Correlation in Time 

Let us now look at correlations in time (Fig. [6|). Al- 
though superconductivity docs not affect the depth of 
the antibunching dip and the lateral coherence as is clear 
from Fig. [4l it influences the correlation function in time 
even at 6* - [Fig. EKa)]. The analytical formulas (|4:27l) 
and (|4.29p show that the correlation function Q at 6* = 
decays as a function of the delay time t ^ t2 — ti like a 
Lorentzian ~ (1 + E^t^)~^ for a normal emitter, while 
formula (|4.30|) reveals that the decay for the supercon- 
ducting emitter is like ^ (Ar)~^ with a different power 
tail [sec (jPlOp for the asymptotic behavior of Ki,{z)]. 
Both cases are compared in Fig. [71 This prolonged tail is 
due to the divergence in the spectrum at the edge of the 
gap (Fig. and is an effect of the superconductor. 

At TT, where the superconductivity of the emit- 
ter prominently manifests itself, the correlation function 
exhibits oscillatory decay in time [Fig.jGjb)]. This is gov- 
erned by the convolution in (|4.3ip : the damped oscilla- 
tion e*'*+'^'/^ in (|4.33p is convoluted with the slow decay 
of |iJ^^^(|Ar|)| - |Ar|-i/2 in KTBi [see ^ for the 

asymptotic behavior of Hi {x)]. The time scales of the 
damped oscillation are given by (|4.25ap . 



iK+/2 = n 



^ijkpr — 2v'uP' jr"^ 
1 + {^2vkFw'^lrf' 



(5.2) 



See Fig. [51 which demonstrates how the parameters w 
and r control the damped oscillation. 



VI. ENTANGLEMENT 

We now turn our attention to the state corresponding 
to the antibunching and bunching correlations. As al- 



ready mentioned in Sec. IIV Al a pair of electrons emitted 
through the Andreev process forms a singlet spin state, 
reflecting the singlet state of the Cooper pair. It is in- 
teresting to explore the possibility that a superconductor 
be a source of entanglement i^iii^i^iiiSiiSiiiS 

Let the two angles (j], Q specify the orientation of 
the spin of an electron, the polar and azimuthal angles 
with respect to a certain quantization axis of the spin. 
The probability that an electron is found at (ri, ti) with 
its spin oriented in the direction (771 , Ci ) and another at 
{ri,ti) {t2 > h) with (??2,C2) is proportional to 

Pr„c.,mCi(2;l) = (<o(l)^i2C2(2)V',BC2(2)^r,iCi(l)>> 

(6.1) 

where 



71 

V'f (r , t) cos — + ^i{r, t)e~^^ sin 



(6.2) 



is the field that annihilates an electron with its spin ori- 
ented along (77, C,). This two-particle distribution function 
is cast into the following form through Wick's theorem: 

PmC2,'niCi (2i 1) 
= 7(2; 2)7(1; 1) 



l7(2;l)|^ 
lx(2;l)|^ 



cos — COS — 

2 2 



-C2) ■ m ■ m 

sm — sm — 



. vi V2 'ni ■ ^2 

sm — cos e^' cos — sm — 

2 2 2 2 



(6.3) 



All the matrix elements of the density operator g of the 
spins of a pair of electrons found at (ri,ii) and {r2,t2) 
(^2 > ii) are reconstructed by combining Pj72C2.r;iCi (2; 1) 
with appropriate sets of the angles to yield 



^(x (7(2; 2)7(1; l)-|7(2;l)|'j 1 

2(|7(2;1)|2-H|x(2;l)|2) \^-){^- 



(6.4) 



where l^-") = (|ti) - |it))/V2 is the singlet state. Note 
that the normalization factor is given by the two-particle 
distribution p^'^\2] 1) in (g^. 

The second term of (|6.4p is the singlet state, while the 
first one is the background that blurs the entanglement. 
At ^? = (when the two detectors are at the same place) , 
the coincident detections with ti = t2 eliminate the back- 
ground and a pure singlet state is extracted. This is be- 
cause the antisymmetric wave function (the triplet con- 
tribution) docs not trigger two detectors located at the 
same point, where the antisymmetric wave function is 
vanishing. In this way, the triplet components arc ruled 
out by the coincident detections and the singlet compo- 
nent is extracted. The entanglement for ~ is just due 
to the symmetry of the wave function and is not supplied 
by the superconductor. 

A remarkable feature of this entanglement is that it re- 
quires no interaction between the electrons. See Ref. [l^ 
for a similar mechanism to extract entanglement, mak- 
ing use of the Bose/Fcrmi symmetry of two-particle wave 
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A| = 0.003 [fj] 




(a) 
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T =t2-tl [h/l^l] 

FIG. 7: (Color online) Normalized coincidence Q vs delay 
time T = t2 — ti for 9 = 0, for normal and superconducting 
emitters. The parameters are the same as those in Fig. |4l 



function. For electronic systems, the interaction-free 
sources of electron-hole entanglement in solid are pro- 
posed in Ref. [H (although the mechanism is different 
from the present one). 

Let us look at the concurrcnco^^ of the state g in (|6.4[) 
extracted by the two detectors. It is given by 



V '27(2; 2)7(1; l)-|7(2;l)P + |x(2;l)P^ 

(6.5) 

This actually gives C{g) = 1 at 6* = 0. It is possible to 
show that the concurrence (|6.5p is nonzero, C{g) > 0, if 



g < 3/4, 3/2 < Q 



(6.6) 
and 7 is 



[note that x is essentially vanishing for 
essentially vanishing for 6 it]. 

The electron pair can even violate Bell's inequalityi^ 
A nonlocal correlation between the two electrons is 
disclosed if the Clauscr-Hornc-Shimony-Holt (CHSH) 
inequality^ 
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(6.7) 



500 1000 1500 2000 

T = t2-tl [h/l^l] 

FIG. 8: (Color online) Normalized coincidence Q vs delay 
time T = t2 — ti for = n, (a) with different values of w and 
(b) with different values of r. Other parameters are the same 
as those in Fig.Hwith |A| = 0.003 [fi]. 



is violated with a certain set of angles t]i^2 and ?][ 2, where 
S = E{r^,,rj2) + E{in,V2) - E{v[,V2) + E{v[,V2) (6-8) 
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and 



P')l.'72 "I" P»?l+7r,')2+7r Prii.ri2+-n P?)i+7i 



Pr]l.7]2 ~^ Pr]i+7T.r]2+7T P7]i,r]2+7r Pr]l+7T.7]2 

(6.9) 

with Prii,rj2 ~ P-qiC-mC- "^'^ present case, the two- 
particle distribution function PriiCi.ri2C,2 (|6.3p yields 



E{r]i,r]2) = -Dcos{iji - 772) 



with a factor 



|7(2;1)P + |X(2;1)P 
27(2; 2)7(1; 1)-|7(2;1)P + |X(2;1)P 



(6.10) 



(6.11) 



Notice here that £^(771, 772) = — cos(7;i — 7/2) represents the 
singlet spin state, and it maximally violates the CHSH 
inequality (|6.7p with e.g. 771 = 0, 772 = — 7r/4, 77^ = 7r/2, 
and 772 = 7r/4. yielding S = — 2V2. The prefactor D {< 
1) in (|6.1ip shrinks the magnitude of S and reduces the 
chance of the violation of the CHSH inequality (|6.7p . The 
CHSH inequality (|6.7[) can thus be violated only when 
D > 1/ \/2. This is accomplished when 



< ^/2/(V2 + 1), \/2/(V2 - 1) < 



(6.12) 



In particular, at 6* = (when the two detectors are at the 
same place), D ~ 1 and the coincident detections extract 
the singlet state, in accordance with the above argument. 

The electron pair for 9 ^ n, corresponding to the 
bunching peak, can also be entangled and can violate 
Bell's inequality, provided the bunching is strong enough. 
In this case, the paired electrons are emitted in opposite 
directions, and this is a more suitable situation for dis- 
cussing nonlocality. This channel is also preferable as a 
source of entanglement, since one need not argue how to 
separate the entangled electrons;^ This is a remarkable 
contrast to the transport setups, in which it is generally 
hard to split the entangled pair. 

The bunching peak at = tt is given by formula (|5.ip . 
There are essentially four independent parameters that 
control the bunching peak, and hence the entanglement, 
i.e. If, r. A, and Eq- See Fig. [21 where the peak and the 
concurrence are plotted as functions of these parameters. 

It is obvious that the gap parameter A plays a signifi- 
cant role for entanglement by enhancing the contribution 
of the Cooper pairs. By increasing |A|, the gap becomes 
wider, and the Andreev emission becomes dominant. The 
enhancement of the entanglement with a wider gap is an 
evidence for the importance of the Andreev emission to 
entanglement. The parameter Ec also works like a fil- 
ter for the Andreev emission. By decreasing Ec, single- 
particle emissions from the quasiparticlc excitations are 
suppressed, and in this way, the background is reduced. 
As a result, the pair emissions become dominant and the 
entanglement is enhanced. 

It is interesting to observe that the size of the emitting 
region, w, also affects the entanglement. The electron 
pairs bunch and are entangled only when they come from 
a small region. This is because the emission from a small 



region ensures that a pair comes from a common Cooper 
pair in the emitter. (If the pair electrons come from 
different Cooper pairs, they are not correlated and do 
not bunch.) That is why the ratio w/^, between the size 
of the emitting region w and the extension of the Cooper 
pair ^, enters formula (|5.ip for the peak of the bunching 
and rules the entanglement of the emitted pair. 

This is essentially equivalent to the effect due to the 
crossed Andreev reflection, which has been extensively 
studied in transport problems ^lii^i^ili^ i i0ii2,24,25 ^j^gj-g 
an electron injected from one normal lead to a super- 
conductor is reflected to the other normal lead as a hole, 
and vice versa. It is shown that the crossed Andreev 
reflections are suppressed as the distance between the 
contacts of the normal leads to the superconductor is 
extended i^i^iii^i^^ It is instructive to rewrite the pair cor- 
relation function xo given in (j4.4p in the following way. 
At far places from the emitter, the first contribution in 
(|4.4p survives and only energies close to the Fermi level 
contribute to xo- By putting Sp- ~ in Sp- + uj^ — 70+ 
(i = 1,2) in the denominators and by substituting the 
integral expression for the tunneling matrix element Tpf, 
given in (j2.7p . the pair correlation xa is cast in the form 

Xa{'ri,ti;r2,t2) 

j d^r[ (fr'^ $(ri ~ r[,h-r2 - r'^,^) 



X 9{r[)g{r'2)T{r[ - r^) 



(6.13) 



with 

$(ri,ii;r2,t2) 



h{Pi)KP2) 



Sp^ -\- Sp2 tO^ 



g-iEpj (tl-t2) 

(6.14) 



and 



^^^^ ^ I (^ 



mkp smkpr 



uJk (27r)"^ kpr 

(6.15) 

for |A| ^ /i. The function !F{r) in (|6.15p is exactly the 
one often found in the arguments on the crossed Andreev 
reflections in the literature f^i^iiiiS'^i^ which describes 
the correlation between the electrons emitted at different 
points r'l and r'2 in the superconductor and decays as the 
emitting points r[ and r'^ are separated far away. The 
emitted electrons are propagated by $(ri — r[,ti;r2 — 
r2,t2) in vacuum, from r[ to ri and to r2, and such 
processes are integrated over the emitting points r[ and 
r2 in the source with a weight function g{r), chosen to 
be Gaussian of size w in (|2.8p in the present analysis. 
This is the physical structure of xo- The formula (|6.13p 
clarifies the connection between the effect of the size of 
the source, w, and the crossed Andreev effects. 

The bunching peak decays as a function of r. It decays 
like ^ {f /kp^^)^^ with oscillations around it [see (|5.ip 
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FIG. 9: (Color online) (a)-(d) The peak of the normalized coincidence Q at 9 = tt and ti = t2, based on the analytical formula 
(|5.ip . The line at Q = 1.5 indicates the threshold for nonzero concurrence C{g) > and the one at Q = V2/{V2 - 1) ~ 3.41 
is that for the violation of Bell's inequality, (e)-(f) The concurrence C'{g) corresponding to (c)-(d). The line at C{g) = 
{3V2 — 2)/4 ~ 0.56 corrseponds to the threshold for the violation of Bell's inequality. The parameters are the same as those in 
Fig. |4] with |A| = 0.003 [n]. The correlation length of the superconductor is ^ ~ 33.8 [\f] for this gap parameter. 



and the asymptotic behavior of Hi, (iz) in (|D9P ]. This 
oscillation originates from the divergences at the edges of 
the gap, E = ±|A|. The reason why the bunching peak 
decays as a function of r is the following. The wave pack- 
ets of the emitted electrons spread as they propagate. 
Even if two electrons are detected at the same distance 
in opposite directions, this docs not ensure that the two 
electrons originate from a common Cooper pair. There 
is an ambiguity to the extent of the spreads of the wave 
packets. The propagator in the free space explains that 
the uncertainty develops up to < y^t/m after an elapsed 
time t, where t is translated into r via t ~ mr/kp- This 



uncertainty should be smaller than ^ for the two elec- 
trons to bunch. This is the reason why the bunching 
peak (|5.ip decays for r ^ kp£^'^. Note, however, that 
this length scale kp^^ is much longer than ^, and the de- 
cay is slow due to the divergence in the quasiparticle 
spectrum. 

VII. ROBUSTNESS 

Let us discuss the robustness of the positive correlation 
at = TT in a non-ideal situation. We consider three types 
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of imperfections: static fluctuations of the diameter w and of the position of the emitting tip, 5r, as well as surface 
roughness that alters the direction of emitting electrons. The shift of the emitting center can be accounted for by 
replacing the tunneling matrix elements in p.5p by 



while the surface roughness would be represented by the deviation of the angle 6 from tt, i.e. 
positive peak at = tt is given, instead of (|5.ip . by (Appendix [E| 



(7.1) 

TT -6. Then, the 



i2Kim/Ec) 
X e 



dw / cr{5r) / d{69)P{w,5r,Se) 



4A/7r 



(7.2) 



where P(w^br^b&) is a probability distribution function 
that characterizes the fluctuations and 8r\\ ^ ri ■ Sr is 
the shift of the emitting center parallel to the directions 
to the detectors fi = —f2- If the emitting center shifts 
toward one of the two detectors, coincident detections at 
the same distance do not imply that the two detected 
electrons originate from the same point, from a common 
Cooper pair, and the correlation is reduced when the 
shift is larger than the size of a Cooper pair. That is 
why iJj'ii/^ appears in the formula. 

When w ^ ^ (which is one of the requirements for 
a strong correlation), the first of the two contributions 
in the absolute value is dominant over the other, as 
mentioned below (|5.ip . Notice that the fiuctuations in 
w and Sr appear there through the ratios 5w'^ /^^ and 
((5r||)^/^^. The positive correlation is hence affected by 
these fluctuations only when they become of order of ^; 
otherwise, it is robust against them. The surface rough- 
ness on the other hand shrinks the height of the posi- 
tive peak mainly by the exponential factor e~''F^ ^^^"^ 
when S0 > 1/kpw. The positive correlation is there- 
fore tolerant of the surface roughness up to an order of 
1/Sk = l/kpSe ~ w. 



VIII. SUMMARY 

We have fully analyzed the correlations, both in space 
and time, of the electrons field-emitted from a supercon- 
ducting tip into vacuum. The superconductivity of the 
emitter leads to positive correlation between the elec- 
trons emitted in opposite directions 9 ^ tt. They can be 
entangled and eventually violate Bell's inequality. A co- 
incidence experiment can directly capture these features. 
Notice that, in contrast to transport setups, one need 
not argue how to separate the entangled pairj^ This is a 
preferable feature, as it makes nonlocality tests easier and 
shows that superconducting nanotips are good sources 



of entanglement. Furthermore, these electron pairs are 
available in vacuum. 

The Andreev emission plays a crucial role. Our ana- 
lytical formulas show that positive correlations and non- 
locality are controlled by the parameters w^/^^, r/kpS,'^, 
and A/Ec- The conditions w^/^^ < 1 and r/kpS,'^ < 1 
ensure that the paired electrons originate from a com- 
mon Cooper pair, and a larger A/i?c makes the contri- 
bution of the Andreev emission more dominant, resulting 
in enhancement of the correlations. Additional require- 
ments, implicitly suggested and/or implied by our anal- 
ysis, would help to enhance these effects and make these 
correlations more manifest and easily observable. For in- 
stance, since the Andreev emission is the only process 
that has a nonzero emission spectrum in the "gap," a 
better correlation would be extracted by selecting ener- 
gies close to the Fermi level. 

Even if entanglement were not observed, the detection 
of the positive correlation would still be a very challeng- 
ing task. This peculiar correlation is a direct manifesta- 
tion of the singlet spin state of a Cooper pair and would 
be a nice probe of the symmetry of the electron pairs in 
a superconductor: if applied to a triplet superconductor, 
for instance, one should find a negative correlation. The 
present setup, field emission into vacuum, would also be 
useful for detecting the anisotropy of unconventional su- 
perconductors. Extension of the present analysis to more 
general superconductors is an interesting future subject. 
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APPENDIX A: GREEN'S FUNCTION UP TO 
THE SECOND ORDER 



Let us carry out the inverse Laplace transformation 
TS]) to obtain Gpp'{t) up to second order in A. 
First, we need to invert the "matrix" in (|3.6p . which is 
formally done as 

g{s) = [P(s)+A2/C(s)]-i =p-i(s)[l+A2/C(s)p-i(s)]-i, 

(Al) 

where 'Dppr{s) = (s + i£p)5^{p — p'). Notice here that it 
is possible to show that 



Det[l + A2/C(s)p-Hs)] 7^ 



(A2) 



for a sufficiently small Ai^ Indeed, the determinant is 
evaluated as 

Det[l + \^lC{s)'b-\s)] = e^'°s[i+^''^(«)^"(^)l, (A3) 
which reads, in the weak-coupling regime. 



l + A2Tr[/C(s)p-i(s)] +0(A4) 



= l + A^y cPpiY[K.pp{s) 
^l + X^ [ d^p [ d^k 



1 



s + 

\TpkUk\'^ 



+ y / d^p / d^fe 



(s + iep){s + iuk) 
^ \TpkUk\^ 
(s - iep){s - iujk) 

\TpkVk\'' 
[s + iep){s - iuJk) 

I \TpkVk? 
[s - i£p){s -f iuJk) 



(A4) 



For a regular tunneling matrix Tpk [i.e., if the spectral 
functions J„(£, w) = / d^p J d^k \TpkUk\^5{ep — e)S{uJk — 
uj) and Jt,(e,w) have good spectral properties: free from 
a nonlocal spectrum like d{e—uj) and vanishing for e, w — *■ 
— ^, oo], the integral in (|A4p is bounded for any s (even 
on the imaginary axis of the complex s plane), and the 
determinant (|A4p can always be made non-vanishing by 
choosing a sufficiently small A. It is therefore possible to 
expand the second factor in the right-hand side of (jAip 
in a power series in A^, 

Gis) = ■D-\s)[l - A2/C(s)p-i(s) + OiX% (A5) 



Gpp' (s) = 



1 



s + i£. 



-5\p-p') 



s + i£p' 



+ 0(A4). (A6) 



Its inverse Laplace transformation yields Gpp' (t) up to 
second order in A, 



Gpp'it)^e~''''''5^ip-p' 



ds 
2Td 



{s+iep){s+iep,) (s+i£p)(s-iep,) , 

e 



{s — iep}{s+iej,i) (s — iep)(s — iSp/ ) / 



0{X% 
(A7) 



Let us take the stationary limit t ^ oo in the inter- 
action picture Gpp'{t) = Gpp'{t)e^^p'*° . It proceeds, for 
instance, as follows: by noting 



1 



{s + iep){s - iSp') 
1 



1 



1 



i{sp + epi)±0+ \s — lEp' s + iEp 
one gets, for t, — > cx) keeping t — to finite, 

ds kii,{s) 



(A8) 



Cb 



i-Ki (s + iep){s — iSpi) 



Cb 



1 



=.(«-»ep')t 



i{ep + £pi) + 0+ s - lEp 



i(£p + Epi) + 0+ S + i£p 



gi£p/(t-to) 



kll,{iep.+Q+)^ 



e p 



(t-to) 



i{£p + £p>) + Q+ 

Recall the formula 



lim 



(A9) 



(AlO) 



Similar treatments are applied to the other components, 
and one ends up with 



APPENDIX B: FAR-FIELD SPECTRA 

Let us sketch the derivations of the energy spectra in 
the far-field limit. Let us first look at the asymptotic 
behavior of the wave functions in (|4.5p .— At far places 
from the emitter, kpr ^ 1. only momenta oriented along 
±fi (z = 1,2) contribute to the propagation of the emit- 
ted electrons. The saddle-point approximation (method 
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of steepest descent) for the integrations over the orienta- where p± are defined in (|4.7p and formula (jA10[) is used, 

tions of the momenta yields One therefore gets (|4.6p . By plugging these expressions 

together with the tunneling matrix (j2.7p and (|2.8p into 

d P Tpk ^ip.r the correlation functions (|4.ip and (|4.3p . and by per- 

y^(27r)^ i £p T i^/c — *0+ forming the Fourier transformation with respect to time 

-, „oo f T ~ T = ti ^ t2i the energy spectra of the correlation hmc- 

/ dppi (P^)fc e'P'- tions in the NESS a nd in the fa r-field limit, ((4A4a|) and 

Jo \epTuJk- «0 (|4l4bl) with and (|4l5b)) . are obtained. 



'271 r 



T{-pf)k -ipr 



dpp- 



27r0(^±Wfc)T( -)fc— — (Bl) Similar treatment to (|Bi p applies to (|4.4|) . yielding (for 



ri = r2 = rj 



Xo{ri,ti\r2,t2) 

/■ ,3, Z""" , /"^ , '^(piri)fcT'(p r2)(-fc) /I 1 



X g'i(Pi+P2)i-g-iepj (ti-t2) 

\2 /• poo poo rp .,,71, . ,,pi(pi+P2)r 

^ ' d'kukvj dp.pj dp2P2 , ^(e— Pi(*i-*^)_e--'=(*i-*=)). (B2) 



In the second term, the poles of the denominators (on-shell contributions) are significant for large kpr [see formula 
(jAlOp ] but such contributions are suppressed by the last factor. The first term is therefore responsible for the far 
fields of the emitted couple. The pole of its first denominator (which describes the propagation of the emitted pair in 
the vacuum) yields 

Xo(r-i,ti;r-2,i2) I d^kukVk H df^p-fc^T(p,,)fcr(p,,)(_fc)e^(P+^^)'^e— ''(*^-*^) + (1 ^ 2), (B3) 

where p = p{—ep) with p{E) defined in ()4.7|) . Its Fourier transformation with respect to r = — ^2 gives the spectral 
representation 

Mri,r2;E) = -'t^e{^,^\E\) J d^ku^Vk ^^2 Jjt^ r[p(-g)f,]fer[p(i^)f2](-fc)e'''^^''^+^^"''^^'' 



2^A2m2 



2^A2m2 











1 




k ^_ 


LUk H 






1 





- b(-E)i-2](-fc)'^ 



(iJ)r2lfce-^[P(^)-P(-^)l'-. (B4) 



Only the first term survives for large kpr, due to formula (|A10[) . 



For the tunneling matrix (|2.7p with the spherical setup 
(Hm), this yields (|4.14c|) with (|4.15cp . by noting the for- 
mulas for the error functions, (|D13p and (|D14p . 



APPENDIX C: ANDREEV EMISSION 



We briefly recapitulate the field-theoretical description 
of the Andreev emission process via a wave operator. The 
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wave operator W is the operator that describes the scat- 
tering of an initial eigenstate of Hq into a scattering state 
by a scattering Hamiltonian H. It is formally defined by 

W = lim e~'"*e'"°*, H ^ Hq + V. (CI) 

t — >oo 

Let \E) be the eigenstate oi Hq belonging to its eigenvalue 
E and suppose that the Hamiltonian H does not have any 
bound state. Then, the wave operator W is given by 

I 



Assume further that Hq admits only one discrete eigen- 
value Eq besides a continuous spectrum E (as is often 
the case in field theory; the vacuum state), and 



(^oll^l^o) =0. (C3) 

Equation (jC2p is valid also for such a discrete spectrum 
of Hq, and the matrix element for the scattering of \Eq) 
into \E) is described by 



(EmEo) = ^— — ({E\V\Eo) + {E\V- ±——V\Eo) + {E\V- ^—^V- ±——V\Eo)) . 

Eq — E + tO+ \ Eq — Hq + zO+ Eo — H + iQ+ Eq — Hq + «0+ / 

, m 



Let us apply this formalism to the present case, i.e. 
to the emission of electrons from a superconductor. The 
Hamiltonian (|2.9p meets the above requirements. In par- 
ticular, the only discrete state is the ground state, vac- 
uum outside and the BCS state inside the emitter, 

|0) - |0)v«)|BCS)s, (C5) 

with Eq = 0. We are interested in the emission of an 
electron pair from this ground state, in the absence of 



the quasiparticle excitations in the emitter. Hence, the 
relevant matrix element is 

{P2S2;Pisi\W\0), |P2S2;Pisi) =4,^4,J0), (C6) 

where, in the final state, two electrons arc found outside 
while the emitter remains in the BCS state. The lowest 
nontrivial contribution to this matrix clement appears at 
the second order in Ht given in (|2.5p and reads^i^^ 



(P2S2;P1S1|W|0) = — ^ 7T-(P2S2;PlSl|i?TT^-^i/T|0) + 0(A4) 

£pi + £p2 ^ Hq — «0+ 

(C7) 



Clearly, the emitted pair of electrons is in the singlet 
state, and the amplitude is symmetric under exchange 
between pi and p2. In principle, one can compute any 
higher-order processes according to (jC2p and (jC4p . It 
is interesting to note that virtual processes (propagators 
between Hx^s) are involved in the Andreev process. 



1. An Integral with the Fermi Distribution 
Function 



An integral with the Fermi distribution functioni^ 



dx ■ 



j/^^ + 1 sin(7ra/^) 



(Re/3 > Rca > 0). 



(Dl) 



APPENDIX D: INTEGRAL FORMULAS 



It is useful to quote some integral formulas used in the 
derivation of some analytical results. Some of the equa- 
tions that follow require extensions of known formulas. 



Bessel Functions 



The density of state of t he superco nductor, which di- 
verges at the gap with 1 / ^JE"^ — |Ap, yields Bessel func- 
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tions, 



f^V{y + \l2) Jo 



dtil^f y-^/^ cos zt 

(Rei/>-l/2), (D2a) 



/¥r(l/2- t/) ^ ' 

{x > 0, iRei'l < 1/2), (D2b) 



/^r(j/ + 1/2) Jo 



dt (1 -t'^y-^/^ cosh zt 



K,{z) 



Hl^\iz) - \/— e"'"<'"+'^"/^ (Re^ > 0, \z\ ^ oo), 

(D8) 

2 " 



mz 



z+(2i^+l)iri/4 



6l(Imz)(l + e2'''^'')e-^-(2'^+i)W4^ 

(Re z > 0, Ini z — > ±oo), 
(D9) 

(Rcz > 0, \z\ -> oo), 



ir(:.)(|)"' (^>o, X-0+). 



(DIO) 



(Rez^> 1/2), (D2c) Note also the fundamental relation 



poo 



K,{z) 



2'^r(i. + l/2) A 

(Re 2 > 0, Rei^ > -1/2). (D2d) 



In particular, the following formulas are used in the text: 
(i) for a real number x. 



0+iO+ 



eosx^^ 7r^(2) 



where 



iHl'\z) = Mz)+iYUz), 
\hI''\z) = J,{z) ~ iY,{z) 

are Hankel functions; (ii) for Rez > 0, 
dt- 



(D4) 



0+ 

du ■ 



2 Jo ^u(l — u) 2 



\Ju{u " 1) 



2 i_i 2 7i 



==ie-^/2(i^o(z/2)-7rz/o(^/2)) 
= ^e"^/^i?^^'(iz/2) (Rez>0), 

by noting the relations 



(D5) 



h{z) = e-'^^'^'^Miz), K,{z) = |e(-+i)-/2i/(i)(z2). 

(D6) 

The asymptotic behavior of some of the Bessel func- 
tions is useful for the discussion in the text: 



V TTX 



x" —x-'' K^{x) = -K^+i{x). (Dll) 



3. Error Functions 

The complementary error function erfc(z) is an entire 
function with no branch cut and endowed with a property 

erfc(-2) = 2 - erfc(2). (D12) 

For Re z > 0, the following expression is available: 



erfc(2) = ±- 



dt- 



t^iz 



(Rez>0). (D13) 



The complementary error function is related to the imag- 
inary error function by 

erfc(z) = 1 + ierfi(iz), erfL(— z) = — crfi(z), (D14) 

which asymptotically behaves as 



erfi(z) ^ le(Imz) + 



1 + 0(1 



(kl-oo) 

(D15) 



with the convention e(0) = for the sign function. 



APPENDIX E: CORRELATION FUNCTIONS 
WITH A SHIFTED EMITTING CENTER 

If the emitting center is displaced in space by a vec- 
tor Sr, the transmission matrix Tpk acquires a phase like 
(|7.ip . In the far- field regime kpr ^ 1 with a small dis- 
placement Sr <^ r, this just induces simple phases in the 
spectra F, X^h, and Xq in (|4.14p of the two-point corre- 
lation functions 7, Xth, and xo- 



(D7) 
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and a similar modification for xoi where p{E) is defined in ()4.7p . 
The same treatment as the one for 7 yields 

7(n,ti;r2,t2) = e''=-('^-'^>''-7(n,ti;r2,t2), (E2) 

while Xth/o is estimated to be Xth/o with substitutions r ^ r — (ri +f2)-5r/2 and ti — i2 ^ ti—t2 + kp{ri—f2)-dr/2fi. 
For the coincident detections ii — i2 = with kp{Sr)'^ /r ^ 1, the integrals involved in Xth/o estimated as [defining 
a = kpiri - 7-2) ■ (Sr/2/i] 

dE ^ ,-^E-/2.±^^Ea ^ P ^ 1 - i£2^2 ^ . . ^ g-*ijV2«± 



V/S2_ |A|2 v/^2_ |A|2 V 2 



^Ag-»|A|V4«±g-|A|^aV4^(2)(^|A|2/4^^ + |A|2aV4), (E3) 



and one ends up with 
■nA 

4rsin(6l/2) 



^ir jkpm 



.g-kpiv COS (U/ Z) ^ih:pl(ri-r2)-dr\ /4r (E4:) 



instead of (|4.35p . Its ratio to (jE2| at 6* = tt gives the integra .nd of ()7.2p 
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